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Abstract 

The nucleon spectral function in infinite nuclear matter is calcu- 
lated in a quantum transport theoretical approach. Exploiting the 
known relation between collision rates and correlation functions the 
spectral function is derived self-consistently. By re-inserting the spec- 
tral functions into the collision integrals the description of hard pro- 
cesses from the high-momentum components of wave functions and 
interactions is improved iteratively until convergence is achieved. The 
momentum and energy distributions and the nuclear matter occupa- 
tion probabilities are in very good agreement with the results obtained 
from many-body theory. 

PACS numbers: 21.65.+1, 24.10.Cn 
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A longstanding problem of nuclear many-body theory is the question to what 
extent short-range correlations are contributing to the properties of nuclear 
matter. While the bulk properties of nuclear matter are mainly affected 
by long-range mean-field dynamics, the picture, however, changes if nuclei 
are probed at large energy and momentum transfer. The spectral functions, 
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obtained for example from A(e,e'p)X experiments are in energy and 
momentum much wider spread than predicted by mean-field dynamics. Ob- 
viously, the high momentum and energy structure of spectral functions has 
important consequences for dynamical processes, e.g. sub-threshold particle 
production on nuclei. 

An overall measure of short-range correlations in nuclear matter is the de- 
pletion of occupation probabilities in ground state momentum distributions 
by about 10%. The processes behind this number are such that states from 
inside the Fermi sphere are scattered into high momentum configurations 
which clearly are not of mean-field nature. As a result, a momentum distri- 
bution with a long high momentum tail is generated extending much beyond 
the Fermi surface. An important finding is that the magnitude and the shape 
of the high momentum component is almost independent of the system under 
consideration while the inner parts, especially in light nuclei, are affected by 
the shell structure and finite size effects. Hence, the high momentum tails 
of the spectral functions are likely to reflect a universal property of nuclear 
many-body dynamics at short distances. 

Theoretically, many attempts have been made to understand short-range 
correlations in nuclei. Approaches based on nuclear many-body theory up to 
explanations referring to the QCD aspects of strong interactions @ have been 
proposed. Obviously, before conclusions on non-standard phenomena can be 
drawn the many-body theoretical aspects of short-range correlations must be 
understood in detail. In fact, the results obtained from many-body theory 
are describing the available data rather satisfactorily. In recent years the 
theoretical results have converged, at least for the depletion of nuclear matter 
occupation probabilities. The majority of the model calculations are using 
Brueckner and Dirac-Brueckner techniques, see e.g. 0, || [|, [6|, [/J. In || a 
correlation dynamical treatment was applied. Most of the approaches use the 
quasi-particle approximation, i.e. a sharp energy distribution for the nucleons 
is assumed (see e.g. 0). Occupation probabilities in finite nuclei could be 
well described in a second RPA approach |I| and by polarization self-energies 
H| , respectively. The Dirac-Brueckner calculations in J7J, including hole- hole 



propagation, led to an extended and numerically rather involved energy- 
momentum structure of self-energies. However, the net effect on binding 
energies and occupation probabilities was only moderate, but improving the 
agreement with empirical data. 

Here, we investigate spectral functions in nuclear matter by quantum 



transport theory [12, [L3|, thereby taking up a proposal of Danielewicz and 
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Bertsch jTJJ . Indeed, the present study is motivated by our recent implemen- 
tation of off-shell effects in a transport theoretical treatment of heavy-ion and 



other nuclear collisions [p~5] , |16|| . A sound theoretical basis how to treat off- 
shell effects in transport equations has been given in |L7| (see also Jf8|). A 
central result of transport theory is that collision rates and correlation func- 
tions are directly related: The calculation of either of the two quantities 
requires the knowledge of the other one. Theoretically, this corresponds to 
a rather involved self-consistency problem for which a direct solution ap- 
parently does not exist. Similar to 0, a practical approach is obtained by 
an iterative procedure. Successive approximations for self-energies, spectral 
functions and collision integrals are obtained be re-inserting the correspond- 
ing quantities from previous cycles of the calculation until convergence is 
achieved. The method is discussed below. Results are presented and com- 



pared to the work of Benhar et al. [fEl who calculated the nuclear matter 



spectral function within the framework of correlation-basis theory and of Ciofi 



degli Atti et al. [20, 21 1 who derived a global parameterization of spectral 



functions for finite nuclei and nuclear matter. 

In quantum transport theory [|l^, [13] dynamical processes are described 
by the one-particle correlation functions 

<7>(1,1') = -z(*(l)*t(l')) 

= z<*t ( r> ( i)), (i) 



g K 



where \l/ are the nucleon field operators in Heisenberg representation. They 
account for the non-stationary processes which introduce a coupling between 
causal and anti-causal single particle propagation. In other words, states 
from below and above the Fermi surface are dynamically mixed as discussed 
above. Correspondingly, in an interacting quantum system the single parti- 
cle self-energy operator includes correlation self-energies S <> which couple 
particle and hole degrees of freedom [|T^, |T3[. Obviously, g K> and S <;> are 
closely related to each other and must be determined by those parts of the 
fundamental interactions producing non-stationary effects. The wanted re- 
lation is obtained from transport theory. After a Fourier transformation to 



energy-momentum representation one finds for the self energies flT2 



v . . / d 3 p 2 du 2 d 3 p 3 du 3 d 3 p 4 du 4 4 4 



(2tt) 4 (2tt) 4 (2tt) 4 
x g < (u 2 ,p 2 )g > (uj 3 ,p 3 )g > (uj 4 ,p 4 ) (2) 
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S < (o; ) p) = ^ (27r)4 (2?r) 4 (27r)4 (2vr)MP + P 2 -P3- P4 )|.M| 2 

x g > (u 2 ,p 2 )g < (uj3,P3)g < (u 4 ,p i ), (3) 



where g = 4 is the spin-isospin degeneracy factor and \A4 1 2 denotes the square 
of the nucleon-nucleon scattering amplitude, averaged over spin and isospin 
of the incoming nucleons and summed over spin and isospin of the outgoing 
nucleons. Since both g <y and S <> describe the correlation dynamics, the 
spectral function can be obtained from either of the two quantities as the 
difference over the cut along the energy real axis. In terms of the correlation 
propagators, 

a(u,p) = i (g > (uj,p) - g < (uj,p)) . (4) 

In a non-relativistic formulation the single particle spectral function is 
found explicitly as 

a \ UJ ,p) = - 3 ——7 NXO , r— r, (5) 



a; 



2mjv 



including the particle and hole nucleon self-energy S, which accounts for 
long-range mean-field dynamics. The width of the spectral distribution is 
determined by the imaginary part of the self-energy, 

V(u,p) = 2ImE(w,p) = i(E > (a;,p) - E < (u,p)). (6) 

From Eqs. (0), (|j) it is apparent that the high momentum, i.e. short range, 
components of nuclear interactions are of primary importance for the energy- 
momentum spreading of the single particle strength. Obviously, the ^-shaped 
quasi-particle distribution is recovered for ImS — ► 0, i.e. vanishing correla- 
tions. 

The correlation propagators g K> are given by 

g < (u,p) = ia(u,p)f(uj,p), (7) 
g > (u 1 p) = -ia(u,p)(l - f(u,p)) (8) 

in terms of the energy-momentum phase space distribution function f(u,p). 
Since we are dealing with a system at T = 0, / reduces to 

f{u,p) = e{u F -u>) (9) 
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with the Fermi energy top. Therefore, 



£ > (u;,p)=0, T(u,p) = — iS < (u;,p) for u> < uj f (10) 
£ < (u;,p)=0, T(u,p) = iS > (a;,p) for uj>uj f . (11) 

In order to obtain T, we have to calculate the self-energies S > and S < . 
Since these quantities themselves depend on T via the spectral function 
a(uj,p), the calculation has to be done iteratively. 

Diagrammatically, the particle- and hole-type transition rate S > and S < , 
Eqs. (0) and (|3p, are of two-particle-one-hole (2plh) and one-particle-two- 
hole (lp2h) structure, respectively. In this respect, they are of the same 
basic structure as the polarization self-energies considered in many-body the- 
oretical descriptions. However, while in many-body theory the polarization 
self-energies are typically included perturbatively in lowest order only by 
performing the integrations over intermediate 2plh and lp2h states with 
quasi-particle spectral functions, e.g. in Ref. fl9| and also Refs. |10|, |TT 



we 



apply a more extended scheme accounting for higher order effects. 

A first attempt to extend the many-body scheme to higher order calcula- 
tions was made in JJJ. Since, as expected, the computational effort was found 
to increase considerably approximations on spectral functions of particle in- 
termediate states had to be invoked. An important aspect of the present 
transport theoretical approach is the fully self-consistent treatment of parti- 
cle and hole strength function in each stage of the calculation. Within our 
iterative approach this is achieved by calculating the correlation self-energies 
with the self-consistently obtained spectral functions. Hence, the converged 
results include a non-perturbative summation of a whole series of np mh 
diagrams. 

Dynamically, the strength of the transition rates is determined by the in- 
medium (off-shell) nucleon-nucleon scattering amplitude Ai. We now make 
the extreme assumption of a constant transition amplitude, as is appropri- 
ate for the short-range part of the NN interaction. Neglecting the energy- 
momentum dependence of the amplitude clearly simplifies the computations. 
Then, M. contributes only as an overall multiplicative factor to the transition 
rates which, at a given density, may be treated as an empirical parameter 
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M. The self-energies E> and E< are then given by 



\ A -\ M \ 2 [ j f j f j 2 f j 2 dcOS $2 , v 

E (u;,p) = ~ 4i ^p y y du 2 J dp 3 p 3 J dp 2 p 2 - - a{uj 2 ,p 2 ) 



xf(u 2 ,p 2 )a(u; 3 ,p 3 )(l - f{to 3 ,p 3 )) I dp 4 p 4 a(uj 4 ,p 4 )(l - f(uj 4 ,p 4 )) 

(12) 



and 



v^</ \ A -\M\ 2 f, f, f, 2 f, 2 fdC08-& 2 , \ 

E (w,p) = y rf^ 3 y gL> 2 y apsPa y d V2V 2 J - — a{^2,P2) 

x(l - f(v 2 ,p 2 ))a(uj 3 ,p 3 )f(u; 3 ,p 3 ) J dp 4 p 4 a(io 4 ,p 4 )f(io 4 ,p 4 ) 

(13) 

with p tot = \p + p 2 \, where four integrations due to the delta function and 
two over the azimuthal angles tp 2 ,tp3 already have been carried out. The 
remaining six dimensional integrals were calculated on a (c<j,p)-grid with 
boundaries < 0.5 GeV and p < p max = 1.25 GeV and a mesh size of 5 
MeV in both directions. 

In the present study we work in infinite nuclear matter at equilibrium 
with the Fermi momentum pp = 1.33 fm -1 , density p=0.16 fm~ 3 and a bind- 
ing energy per particle of up = —16 MeV. The real part of the self-energy 
ReE is chosen to be independent of momentum and energy. In fact, a con- 
stant ReE only serves to define the scale for the excitation energy uj which 
in our case is given by w > up. The energy and momentum dependent 
pieces of ReE play a more important role because they will modify the pole 
structure of the propagators. Inside the Fermi sphere typically a quadratic 
momentum dependence is found leading to a scaling of the kinetic energy by 
an effective mass and a global compression of spectra will be found. Hence, 
leaving out the momentum dependence will affect the final results only to 
a minor degree. Neglecting the energy dependence might introduce larger 
uncertainties because this amounts to a violation of analyticity as expressed 
by the dispersion relation between real and imaginary parts of polarization 
self-energies. Since we do obtain an energy dependent imaginary part analyt- 
icity could in principle be restored by calculating ReE dispersively. However, 
using a constant, energy independent transition matrix element clearly fails 
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for large positive energies because the imaginary part would continue to in- 
crease beyond any limit with the level density of two-particle states. Thus, 
the imaginary part of the self-energy is not available over the full range of 
energies necessary for solving the dispersion equation. Rather than introduc- 
ing an energy cut-off as an additional free model parameter we decide at the 
present stage to completely neglect the energy dependence of ReS. 

The average amplitude was adjusted such that the spectral functions 

obtained with many-body theoretical methods by Benhar et al. [19| are 

/ \ 1/2 

reproduced leading to ( |.M| 2 ) =207 MeV fm 3 . Relating this value to on- 

shell processes would correspond to a constant cross section of about 20 mb. 
In Fig. [I] the full energy- momentum structure of the resulting nucleon spectral 
function 

P{uj,p) = N a(oo,p) (14) 

is displayed, where N is a normalization constant. The same normalization 
as in the work of Benhar et al. [|19j is used. Our results are in remarkable 



agreement with the calculations of Ciofi degli Atti et al. as seen by 

comparing to Fig. 10 in that reference. 

For a more quantitative comparison cuts for several momenta are shown 
in Fig. as a function of energy uj below the Fermi surface and compared to 



the results of Ref . |T9[ . The agreement is surprisingly good considering the 
seemingly drastic approximation invoked for M. and the neglect of analyticity. 

We conclude from this agreement that the spectral function of nuclear 
matter in the hole-sector is rather insensitive to the specific energy and mo- 
mentum structure of the interaction once a fully self-consistent calculation 
is performed. The effect of self-consistency is illustrated in Fig. 2 where 
spectral functions obtained after the first iteration are compared to the fully 
converged results. The iteration scheme was initialized by choosing spec- 
tral functions with a constant width of 0.1 MeV as start value. On a level 
of about 30% global features, as e.g. the location and height of the peak 
structures, are described already by the first iteration. At first sight, this 
seems to support the pertubative approach usually applied in many-body 
theoretical approaches. However, during the self-consistency cycle strength 
is re-distributed from the peak into the tail regions. Hence, higher order 
polarization effects act in a similar way as increasing the effective coupling 
strength. The results show that the full transport theoretical calculation 
apparently leads to a self-consistent adjustment of spectral functions which 
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ultimately is dominated by phase space effects. 



The agreement with the many-body calculations of [19] extends also to 



the single nucleon momentum distribution in nuclear matter defined as 

/U) F 
duP(u,p). (15) 
-oo 

Results are displayed in Fig. |3|. Again the agreement is close to perfect 
except for the increase of the distribution towards pp, which is in contrast to 



the calculations of [ I9f . Calculations in a schematic model using an empirical 



form for T(uj,p) with parameters given in Ref. J23| show that this behaviour 



is a direct consequence of the missing analyticity in the self-energies. The 
schematic approach also confirms that the momentum dependence of ReS 
indeed leads to minor effects in the bulk behaviour of the spectral functions 
except close to the Fermi surface. The momentum distribution is affected 
mainly in the extreme momentum tail where a steeper slope is found. In this 
context it is important to note that n{p) must decrease stronger than C(^r) 
in order to have a convergent result for the kinetic energy distribution of the 
correlated ground state. The transport theoretical results in Fig. |3| fulfill 
this constraint. 

Finally, we remark that the shape of the momentum distribution is closely 
related to the value of Ai: Increasing the value - corresponding to a stronger 
interaction amongst the nucleons - would increase the occupation of states 
above pp and soften the Fermi edge, leading to a better description of the 
data in Fig. |^ for p > p F . 

Summarizing, correlations in nuclear matter have been described by an 
approach that allows to go beyond the perturbative level to which conven- 
tional many-body methods are constrained in practice. In a first application, 
the transition rates were calculated with a momentum independent nucleon- 
nucleon scattering amplitude. That a transport theoretical description with 
a simple approximation for the scattering amplitude could reproduce the 
momentum distribution obtained in state-of-the-art many-body calculations 



had already been observed by the authors of [14]. The results here show 
in addition that also the spectral functions are dominated by phase space 
effects rather than by the off-shell momentum structure of interactions. The 
excellent agreement of the nucleon spectral function at far-off-shell energies 
and momenta with those obtained in sophisticated state-of-the-art nuclear 
many-body calculations also provides some a-posteriori justification for the 
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method to include off-shell effects in nuclear transport calculations developed 
in Ref. pfl. 
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Figure 1: Energy and momentum dependence of the nucleon spectral function 
for energies below lo f . 
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Figure 2: Nucleon spectral function for different constant momenta for ener- 
gies below up- The dotted curves have been obtained after the first iteration 
step using a constant initial width of 0.1 MeV (see text). The dashed curves 
show the results of Benhar et al. [|19l [22]]. 
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Figure 3: Momentum distribution in nuclear matter calculated from the 



spectral function. The dashed curve is the result from Benhar et al. |T9 
The data are from [EI 
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